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Abstract 

The present work devoted to the finding explicit solution of a boundary problem 
with the Dirichlet-Neumann condition for elliptic equation with singular coef- 
ficients in a quarter of ball. For this aim the method of Green's function have 
been used. Since, found Green's function contains a hypergeometric function of 
Appell, we had to deal with decomposition formulas, formulas of differentiation 
and some adjacent relations for this hypergeometric function in order to get 
explicit solution of the formulated problem. 
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1. Introduction 

Partial differential equations with singularities (PDEwS) have numerous ap- 
plications in real life processes [1,2]. Most famous equation of this kind is 
Chaplygin's equation [3], which can be written as 

K{y)u xx + u yy = 

or in some particular cases as 
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Latter one called as Tricomi's equation and was studied well by many authors 
[4-7]. Due to possible applications and natural mathematical interests for gener- 
alization, (PDEwS) were studied intensively and investigations are still going on. 
One can find bibliographic information on this in the monographs [8-10]. Omit- 
ting huge amount works, devoted to investigations of boundary problems and 
potentials for two dimensional cases of aforementioned equations, note works by 
A.V.Bitsadze [11], A.M.Nakhushcv [12], M.S.Salakhitdinov and B.Islomov [13], 
where three dimensional mixed type equations with singularities were investi- 
gated reducing them into two dimensional case using Fourier transformation. 

Dirichlct and Dirichlet-Neumann problems for elliptic equation with one 
singular coefficient in some part of ball were investigated by C.Agostinelli [14], 
M.N.Olcvskii [15]. Recently, I.T.Nazipov published a paper devoted to the in- 
vestigation of the Tricomi problem in a mixed domain consisting of hemisphere 
and cone [16]. We also note works [17-19], where fundamental solutions and 
boundary problems for elliptic equations with three singular coefficient were 
subject of investigation. By J.J.Nieto and E.T.Karimov [20] the Dirichlet prob- 
lem for an equation 

2a 2/3 

Ha,p(u) = u xx + u yy + u zz -\ u x H u y = 0, 0<2a,2/?<l (1) 

x y 

was studied in some part of ball. 

In the present work handling the method of Green's function we find ex- 
plicit solution of a boundary problem with the Dirichlet-Neumann condition for 
elliptic equation with two singular coefficients in a quarter of a ball. 

2. Main result 

We consider Eq.(l) in a domain 

n = {{x,y,z) : x 2 + y 2 + z 2 < R 2 , x > 0,y > 0, -R < z < R} 

which is a quarter of a ball. 

Problem D-N. Find a function u(x,y,z) <G C (f2) flC 2 (!l), satisfying 
Eq. (1) in ft and conditions 

u (x, y, z) \ X=Q = n {y, z) , (y, z) <G fii, (2) 

y 2fi u (x, y, z) \ y=0 = v 2 (x, z) , (i,z)efl 2 , (3) 
u(x,y,z) = ip(x,y,z) , (x,y,z) eE. (4) 

Here 

Oi = Ux,y,z) : y 2 + z 2 < R 2 , x = 0,y > 0, -R < z < R\ , 
£1 2 = {(x,y,z) : x 2 + z 2 < R 2 , x > 0,y = 0, -i? < z < R} , 
§ = {(x, y, z) : x 2 + y 2 + z 2 = R 2 , x > 0, y > 0, -R < z < R} , 

n (y, z) , V2 (x, z) , ip (x, y, z) are given continuous functions fulfilling matching 
condition n (y, z)\ Ti = ip (x, y, z)\ Ti , (Ti := y 2 + z 2 = R 2 ). 
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Figure 1: Domain Q 



Remark. The uniqueness of solution of the problem D-N can be proved 
using classical "a-b-c" method as in the work [20] or using extrcmume principle 
for elliptic equations [8]. We note also that the uniqueness theorem can be 
proved for more general domains. 

First we give a definition of the Green's function of the problem D-N. 

Definition. We call the function G (M, Mq) as Green's function of the 
problem D-N, if it satisfies the following conditions: 

1. it satisfies equation 

G XX + Gyy + G ZZ + —U X + —Uy = S (M , M Q ) 

x y 

in 

2. it satisfies boundary conditions 

GUo = 0, y 2p G y \ y=0 ^0, G|g = 0; 

3. it can be represented as 

G (M, M Q ) = q (M, M ) + q* (M, Mo) • (5) 

Here 5 is the Dirac's delta-function, M(x,y,z) is any and Mo(x,y, z) is fixed 
point of O and Ma (xq, zq, zo) is a point symmetric to Mo regarding the spherical 
surface of considered domain, i.e. 

r>2 r>2 n2 
A / 1 IX 2 2 2 2 

^0 = — T>2 X 0' V0 — — 152^0, z = ~T>2 Z 0> = X + % + Z ' 

it H Q it 

q (M, M ) - k (r 2 )"^~ f (xx ) l - 2a F 2 (^-a + p,l-a,l3;2- 2a, 2/3; f , ^ 

(6) 

is one of the fundamental solutions of Eq.(l) [20], 

Fi (a;bi,b 2 ;c 1 ,C2]x,y) = > a:;/ 

i^o ( c i)«( c 2), i! J ! 
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is Appel's hypergeometric function [21], 

j£ 3-2ct+2,3 



q* (M,M ) 



Rq ) 



q(M,M ) 



is a regular part of the Green's function (5), i.e. satisfies Eq.(l) in any point of 



v 



„2' 



(z - a; ) 2 + (y- y f + {z- z f 



„2 _ („ i „ \2 , /. „. \2 , / „ \2 2 _ / „ \2 , /„. \ 2 i / „ \2 



(7) 



r 1 = (a; + xo) + (y - y ) + (z - z ) , r 2 = (x - x ) + (y + y ) + (z - z ) , 

_ 1 r(l - a)r(/3)r(2 - 2a + 2/3) 
_ 2tt T(2 - 2a)r(2/3)r(l - a + j3) ' 

Excise from the domain Q a small ball with its center at M and with radius 
p > 0, designating the sphere of the excised ball as C p and by O p denote the 
remaining part of Q. 

Applying Green's formula [22] in this domain, one can obtain the following: 



ffx 2a y 2P 

C„ 



dG du 

u— G-^— 

On on 



ds = JJ y 20 n(y, z)G* (M, M ) dydz 

dG W 



+ // x 2a u 2 {x, z)G** (M, Af ) dxdz + jj x 2oL y 2 ^ 'tp(a)— do- 



Here G* (Af, M ) = x 2 



r23 dG(M,M ) 
dx 

outer normal to the dSl. 

Let us first evaluate an integral 



, G**(M,M ) = G(M,M )\ y=0 , n is 



x=0 



dG 



x za y zp u—ds. 
an 



(9) 



Based on (5) we write (9) as follows: 



I = h+h 



2a 2/3 u 1 i 

x y ' u—ds 
on 



2a 2/3 9q* 

x y ^u—-ds, 
on 



where 



dq dq dq dq 

— = — cos(n, x) + — cos(n, y) + — cos(n, z). 



Formal calculation gives 

ox ox 



F 2 (...) + P 1 P 2 



dF 2 | m 

<9£ x dn 



(10) 
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where Pi = (r 2 ) a 13 2 , P2 = (2^0) 1 2a . 

Using formula of differentiation for hypergeometric function F 2 (...) [21] 

d l+J F 2 (a,b 1 ,b 2 ;c 1 ,c 2 ;x,y) 



dx l yi 
("), +J (M, (b2 

(ci)j {02) j 



F 2 (a + i + j, &i + i, b 2 + j; c x + i, c 2 + j; a;, y) 



we get 



3£ 



F 2 ( - - a + /3, 1 - a, /3; 2 - 2a, 2/3; £, 77 
+ /?)(! 



(2 - 2a) 



-F 2 \--a + p,2-a,p;3-2a,2p^,r ] ), 



drj 



a + 0, 1 - a, P; 2 - 2a, 2j3\ £, 77 
{ ~ 2 " ° + '^-F 2 Q- a + /3,l-a,/? + l;2-2a,2/?+l;£,r7 



2/3 



Considering 



<9Pi _ 2 (or - x ) 
<9x r 2 



a-y8- 



SP 2 _l-2a 

ax a; 



4x 2(.t-x ) 2(x-x ) 



we have 




dq 

dx 


kPiP 2 


2(x 


xf 2 ( 


'3 


-PA 


xF 2 ( 




xF 2 ( 


^ -a + /3,2 


xF 2 ( 


5 


-p,2 


xf 2 | 


^ -a + p,l 



1 - 2a 



a + P 



a + P 77 



Now we apply to (11) the following adjacent relation [21] 

x— F 2 (a + l,6i + l,6 2 ;ci + l,c 2 ;a;,y) 

+y — F 2 (a + 1,61,62 + l;ci,c 2 + l;x,y) 
c 2 

= P 2 (a+ l,6i,6 2 ;ci,c 2 ;a;,2/) - P 2 (a,6i,6 2 ;ci,c 2 ;a;,y) 



1 - a 
2 -2a 
P_ 
2P 



(11) 
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and will obtain 

dq _ kP x P 2 
dx r 2 



2(x-x Q )(a-8 



*F 2 ( - - a + 0, 1 - a,0; 2 - 2a, 2/3; £, r, 



+2x 



(a-0- 



+kP 1 P 2 



1 - 2a 



F 2 [--a + 0,2-a,0;3-2a,20;^,r 1 



(12) 



r 2 F 2 [--a + 0,l-a,0;2-2a,20^,r, 



Similarly we get 
dq kP x P 2 



dy 

xF 2 , 
" \2 

+2y (a 



2(y-y )[a-B 



a + 8, l-a,/3;2-2a,2/3;^,ry 



F 2 



(13) 



a + 8, 1 - a, + 1; 2 - 2a, 28 + 1; £, »j 



dq _ kP Y P 2 
dz r 2 



2 (z - zq) (a - - 



-a + 0,l-a,0;2-2a,20;^,r] 



Taking (12)-(14) into account, from (10) we get the following: 
d 



(14) 



dn 



q (M, M ) = [a 



kPiP 2 



xF 2 \--a + 0,l-a,0;2-2a,20;^,r } 



— r lnr 2i , fcp i p 2 

dn r 2 



a-0- 2 



2x F 2 [ --a + 0,2-a,0;S-2a,20;^,r 1 



+2y Q F 2 ^--a + 0,l-a,0 + l;2-2a,20+l;^,r] 

+kP 1 P 2 1 ^F 2 

x 



a + /3,l-a,/3;2-2a,2/3;^,77 . 



(15) 



G 



Now consider the integral 
h= ff x^u^dS = In + I 12 + I 13 

= jj* 2a y 2fi kp ^ 

xF 2 Q-a + /3,l-a,/3;2-2a,2/3;£,^ ^ [lnr 2 ] u (M ) ds 



x 
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2x0^2 ( g -a + /8,2-a,^;3-2a, 2/3^,77 



- 2 M , G ( ^ - a + /3, 1 - a, /? + 1; 2 - 2a, 2/3 + 1; <£, 



u (M ) 

x 2a y 2fi kP x P 2 ^—^F2 Q - a + /3, 1 - a, /3; 2 - 2a, 2/3; £, »?) « (Af ) ds. 

We use the following spherical system of coordinates [22]: 

x = x + psin#cosV', y — yo + psin#sin?/>, z = z a + pcosO, 

<9 <7T, 0<V<27T, < /9 < OO. 



Then we have 





I n = -2 (a - /3 - |) k J dtp J x 1 ^ 2a (x + p sin 9 cos if)) (y + p sin 9 sin V>) 2/3 
o o 

xii (x + p sin # cos?/;, y + p sin # sin?/;, zo + pcos9) (p 2 ) a ^ 1 
xF 2 (| - a + 0, 1 - a, /3; 2 - 2a, 2/3; £ s , r? s ) sin 9d9. 

Let us first evaluate hypergeometric function F 2 . We use decomposition formula 
[21] 



F 2 (a,bi,b 2 ;c 1 ,c 2 ;x,y) 

(ci),^)^! 



E 'tul-u?' 1 ^'^ ( a + *, &i + »; ci + i; x) ■ 2 f\ (a + i,b 2 + i; c 2 +i;y), 



i=0 

after then using auto transformation formula [23] 



2 Fi (a, b, c; x) = (1 - x) 2 F 1 ( c - a, b, c; _ ^ 
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we obtain 



F 2 (^-a + 0,l-a,/3;2-2a,20^ s , Vs ^ = 
~ (|- a + ^(l-a) i(/ 3) 
^ (2 - 2a), (2/3), i! ^ s 

!5 
--a + /? + i,l-a; + i;2-2a + i;£ ; 
5 
--a + /3 + i,/3 + i; 20 + i;r] s 

^ (2-2a) i( 2^i! Cs) (1 ™ 

x 2 f\ ( - \ - a - 0, 1 - a + i;2 - 2a + i; j^— 

x 2 F! ( -^-+a + 0,0 + i;20 + i; ''" 



Here 



2 



= 1 J", ^— = 1 g - , r 2 s = 4a; 2 l + 4x psin#cos tA, 



r 2 — At 2 

t _ I - r -2 ' _ 1 - 2 ' 'is - 

?s 1 r ls 2s 2 2 

HL =4^ + 4?/opsin6isin^, 1-^ = ^,1-^ = ^. 
After some evaluations we find 

F 2 Q - a + 0, 1 - a, 0; 2 - 2a, 2/?; ^ = (r? 8 ) ^ HL)^ (p 2 ) ^ 

2- (2 -2a). (20).*! J , 



x 

i=0 



x 2 Fi I - a- 0,1- a + i;2-2a + i; s 



2 " ' .6-1 

x 2 ^i l~+a + 0,0 + i;20 + i;- ?/ '' 



2 jy s - 1 

Considering [23] 

^ (c)„n! r(c-a)r(c-6) 
c + 0,-1,-2,..., SR(c-a-6) > 0. 

we deduce 

/ 1 p 2 

lim 2-F1 — - — a — 0, 1— a + z; 2 — 2a + i; 1 =- 

p^o V 2 rf s 



r(2-2a + »)r(/3 + |) 
r(|-a + /3 + «)r(i-a) : 
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Now we get 



/ 5 p 2 
\im 2 F 1 i--+a + /3,/3 + i;2/3 + i;l- ^y- 
p^o \ 2 r| s 

T (2/3 + r (-a + I) 

~ r(§-a + /3 + i)r(^)" 



lim I n = -2tt • A; • ( a - /3 - I V 2 2a -^ +1 u (M ) 
p^o \ 2 / 



r(|- a )r(l + /?) 
r(i- a )r(/3) 



where 



" (|_ a + ; 8).(l- a ). (ff). r(2-2a + i)r(2/3 + i) 
45 ^ " (2-2^(2/3)^! " ' r2(|- a + /3 + i ) ' 

Considering [23] 

r (2 - 2a + i) = T (2 - 2a) (2 - 2a) 2 , T (2/3 + i) = T (2/3) (2/3), , T (§) = ^, 
r f5 \ = ^r(2-2a + 2/3)(f-a + / 3) fg_ 

making some calculations we obtain 

2- 2a + 2 ' 3 r(2-2a)r(2/3)r(l-a + /3) 



^ (| - a + /3) T (2 - 2a + 2/3) T (| - a) T (| + /3) ■ 



Therefore 



.. T , n r (2 - 2a) T (2/3) T (1 - a + /3) .... 
iTo 7 " - fc • 2 - r (2 - 2a + 2/3) F (1 — a) T (/3) (M ° } ' 

If we choose fc as in (7), we will have 

lim In = u (M ) . 

p^O 



By similar evaluations one can get that 



lim ii2 = lim 7i3 = lim /2 = 0. 

p^O p^O p^O 



If we consider an integral 

[ [ ^Wr, 

dn 
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using above given algorithm for evaluations, we can prove that 

lim f I x 2a y 2 ^G^dS = 0. 



p^o J J 

c 



dn 



Now we can formulate our result as the following 

Theorem. The Problem D-N has unique solution represented as follows 
u (M ) = JJ y 2/ Vi (y, z) G* (M, M ) dydz + JJ x 2a v 2 (x, z) G** (M, M ) dxdz 



The particular values of the Green's function is given by 

y) x a 



G* (M, M ) = k (1 - 2a) x 1 ^ 2 " 

2 F 1 (l- a +l3.,l3.2l3-r lo:c ) _ ( 2 ~ 4a 2*1 ( f -*+P ,^,2j3;»? 03; ) 

(^U =0 )-i +a -' 3 



(r- 2 |x=o)-2- 



G** (M,M ) = fc(xx ) 1- ' 

2Fi(f-a+/3,l-a,2~2a;g oi ,) 
(^1^0)-*+°^ 



2-4a 



2 Fi(§-a+j3,l-a,2-2a;g 01< ) 
(^1^0)-*+-^ 



Here, 



r 2 | y=0 = (R 

£,0y — — 72J 



4xxq 



yya \ 

R ) 

XXq ^ 

R , 
, VOx 



(R- 

(R- 
iyyo 



zza_\ 
R ) 

ZZQ s 

R , 



y + z o ji 



Oy 



R 2 

AR Z xxq - 



x a+yg „2 _ n2 

^+Vo z 2 -R 2 , 



R 2 r 2 U =0 ' 
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